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Abstract. The aim of this paper is to propose an unambiguous intrinsic formahsm for higher- 
order field theories which avoids the arbitrariness in the generahzation of the conventional de- 
scription of field theories, and implies the existence of different Cartan forms and Legendre trans- 
formations. We propose a differential-geometric setting for the dynamics of a higher-order field 
' theory, based on the Skinner and Rusk formalism for mechanics. This approach incorporates 

aspects of both, the Lagrangian and the Hamiltonian description, since the field equations are 
formulated using the Lagrangian on a higher-order jet bundle and the canonical multisymplectic 
, form on its affine dual. As both of these objects are uniquely defined, the Skinner-Rusk approach 

has the advantage that it does not suffer from the arbitrariness in conventional descriptions. The 
, result is that we obtain a unique and global intrinsic version of the Euler-Lagrange equations for 

' higher-order field theories. Several examples illustrate our construction. 
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1. Introduction 

During the last decades of the past century, there have been different studies and attempts to 
' define in a global and intrinsic way the higher-order calculus of variations in several independent 
0> , variables. The standard geometric picture starts with a Lagrangian function L : J^tt — >■ M 
where J^vr is the /cth-order jet bundle of a fiber bundle tt : E ^ M. The main objectives are 
O I to describe the associated Euler-Lagrange equations for sections of the fiber bundle, to derive 
\Q • Cartan forms for use in intrinsic versions of the above equations, and to construct adequate 
^ . Legendre maps which permit to write the equations in the Hamiltonian side (see, for instance, 
[Ill3IIIl|12lll6l[l9l|25l|271[33l|3l[37]for further information) . 

The situation is well established for the case of one independent variable (higher order mechan- 
ics) and for the case of first order calculus of variations [HI [T7t [T8| I2l] . In this last situation, 
the typical expression of the Cartan form associated in classical mechanics to a Lagrangian 
5^ I L : J^TT — )■ M may be written as S*{dL) + Ldt, where S* is the adjoint of the vertical endomor- 
phism acting on 1-forms. In order to generalize this concept to higher order field theories, one 
needs to define a mapping from 1-forms (the differential of L) to m-forms and to incorporate in a 
global way the higher order derivatives. This is one of the reasons for the degree of arbitrariness 
in the definition of Cartan forms for Lagrangian functions L : J^vr — )■ M, if A; > 1 and dimM > 1. 
In other words, there will be different Cartan forms which carry out the same function in order 
to define an intrinsic formulation of Euler-Lagrange equations. The main reason of this problem 
is the commutativity of repeated partial differentiation. Therefore, the Cartan form is unique if 
(and only if) either or m equals one. 
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In the literature, we find different approaches to fix the Cartan form for higher order field 
theories. A direct attempt is the approach by Aldaya and Azcarraga [H |2]. Another point of 
view is that by Arens [3], which consists of injecting the jet bundle J^ti to an appropriate first- 
order jet bundle by the introduction of a great number of variables into the theory and Lagrange 
multipliers. From a more geometrical point of view, Garcia and Munoz described a method of 
constructing global Poincare-Cartan forms in the higher order calculus of variations in fibered 
spaces by means of a linear connections (see |15[ IT6j). In particular they show that the Cartan 
forms depend on the choice of two connections, a linear connection on the base M and a linear 
connection on the vertical bundle Vt^. Later, Crampin and Saunders [37j proposed the use of 
an operator analogous to the almost tangent structure canonically defined on the tangent bundle 
of a given configuration manifold M for the construction of global Poincare-Cartan forms; this 
operator depends on the chosen volume form on the base. 

In our paper, we propose an alternative way, avoiding the use of additional structures, working 
only with intrinsic objects from both the Lagrangian and Hamiltonian sides. This formalism 
is strongly based on the one developed by Skinner and Rusk [55| HO] . In order to deal 
with singular Lagrangian systems. Skinner and Rusk construct a Hamiltonian system on the 
Whitney sum TQ © T*Q of the tangent and cotangent bundles of the configuration manifold 
Q. The advantage of their approach lies on the fact that the second order condition of the 
dynamics is automatically satisfied. This does not happen in the Lagrangian side of the Gotay 
and Nester formulation, where the second-order condition problem has to be considered after 
the implementation of the constraint algorithm (see [201 ED [22]), besides other formalisms which 
include the second-order condition from the very beginning (see ^8^ ilOj). 

For higher-order field theories, we start with a Lagrangian function defined on J^tt. We consider 
the fibration tiwoM • — ^ ^5 where Wq = J^ti Xjk^i^ A™(J^^^7r) is a fibered product, the 
velocity-momentum space. On Wq we construct a (pre-)multisymplectic form by pulling back the 
canonical multisymplectic form on A^(j'^~^7r), and we define a convenient Hamiltonian from a 
natural canonical pairing and the given Lagrangian function. The solutions of the field equations 
are viewed as integral sections of Ehresmann connections in the fibration Trwo,M '■ Wq — > M. 
In this space we obtain a global, intrinsic and unique expression for a Cartan type equation for 
the Euler-Lagrange equations for higher-order field theories. Additionally, we obtain a resultant 
constraint algorithm. Our scheme is applied to several examples to illustrate our method. 

Apart from the lack of ambiguity inherent in our construction, we emphasize that our formalism 
can be easily extended to the case of higher-order field theories with constraints and optimal 
control problems for partial differential equations. In this way, we obtain a unified, geometric 
description of both types of systems, with possible future applications in the theory of symmetry 
reduction and the construction of numerical methods preserving geometric structure (see |28|). 
This will be the topic of future research. 

While finalizing this paper, we found out about the work of L. Vitagliano jH] who independently 
used the unified formalism framework to study higher-order field theories, using techniques from 
secondary calculus. 

Throughout the paper, lower case Latin (resp. Greek) letters will usually denote indexes that 
range between 1 and m (resp. 1 and n). Capital Latin letters will usually denote multi-indexes 
whose length ranges between and k. In particular, in section ^ and all later sections, I and 
J will usually denote multi-indexes whose length goes from to /c — 1 and to k, respectively; 
and K will denote multi-indexes whose length is equal to k. The Einstein notation for repeated 
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indexes and mult i- indexes is used but, for clarity, in some cases the summation for multi-indexes 
will be indicated. 



2. Jet Bundles 

Let {E, TT, M) be a fiber bundle whose base space M is an orientable differentiable manifold 
of dimension m, and whose fibers have dimension n, thus E is {m + n)-dimensional. Adapted 
coordinated systems will be of the form where (a;*) is a local coordinate system in M 

and (u") denotes fiber coordinates. We fix a volume form rj on the base manifold M. For a 
compatible chart (x*) with respect to the volume form, rj is written d"^x = dx^ A ■ ■ ■ A dx™, and 
we will write d"^~^Xi for the contraction iQ/Q^id"^x {d^~'^Xij = iQ/Q^jd"^~^Xi and so on). 

Given a point p G M, let (p^ip : M — > E be two smooth local sections around p. We say that 
(f) and ip are k-equivalent at p (with A; > 1) if the sections and all their partial derivatives until 
order k coincide at p, that is, if 



d' 

4>{p) = i^ip) cind 



■ ■ ■ dx^'' 



Qk^C 



■ ■ ■ dx'^^ 



for all 1 < a < ra, 1 < < m, 1 < j < /c. Note that this is independent of the chosen 
coordinate system (adapted or not) and, therefore, to be ^-equivalent is an equivalence relation 
(see |[32l EHl [36], for more details). 

Definition 2.1. Let {E,tt,M) be a fiber bundle and p G M. Given a smooth local section 
G rp(7r), the equivalence class of ^-equivalent smooth local sections (with k > 1) around p 
that contains (j) is called the k-jet of (p at p and is denoted The set of all the fc-jets of local 
sections, that is, 

: peM,cPeT,in)}, 

is called the k-th jet manifold of tt and denoted J^vr. 

These sets have interesting structures and relations between them, but before we present them, 
we will introduce a particular multi-index notation. 

Note 2.2 (The multi-index notation, [36j). Given a function / : M"^ — > M, its partial derivatives 
are denoted by 

d'f 



fil 



«2 



Since all the functions that we consider are smooth enough, their crossed derivatives coincide. 
Thus, the order in which the derivatives are taken is not important, but the number of times with 
respect to each variable. 

Another notation to denote partial derivatives is defined through multi-indexes. A multi-index 
/ will be an m-tuple of non- negative integers. The i-th component of / is denoted Addition 
and substraction of multi-indexes are defined componentwise (whenever the result is still a multi- 
index), (J ± J){i) = ± J{i)- The length of / is the sum |J| = X^j-^l^)' factorial 
/! = nj/(i)!. In particular, Ij will be the multi-index that is zero everywhere except at the i-th 
component which is equal to 1. 
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Keeping in mind the above notations, we will denote the partial derivatives of a function 
f -.W — ^ M by: 

^1^1 / ^7(l)+7(2)+-+7(m) J 

Thus, given a multi-index /, denotes the number of times the function is differentiated with 
respect to the i-th component. The former notation should not be confused with the latter one. 
For instance, the third order partial derivative (with / : _^ M) is denoted ^32 and 

/(o,2,i) respectively. 

Let {E,n,M) be a fiber bundle as before. An adapted coordinate system on the total 

space E induces adapted coordinates (with < |/| < k) on the A;-jet manifold J'^vr given 

by: 



dx^ 



from where we deduce that the dimension of J^-k is 



T Tfc , /m - 1 + /\ 

dim J TT = m + 77, ■ > 

1=0 ^ ^ 



It is readily seen that ( J'^vr, VTjt, M) is a fiber bundle, where 

^fc(jp0) = V (ill coordinates Tx^ix!' ^uf 



ix% 



Note that any local section G rp(7r) can be lifted to a local section in TpliTk) defining its lift by 
(see Diagram [1]): 

(jV)(p)=j>. 

It is also clear that the fc-jet manifold J'^tt fibers over the lower order /-jet manifolds J'tt (see 
Diagram H]), with < / < A;, where by convention J^tt = E and where the projections are given 
by: 

'^k,i{ip4>) = jp4> (in coordinates -Kk^i{x\uJ) = {x\uj), with < |/| < A;,0 < \ J\ < I). 
In particular, ( J'^vr, vr^ J'^^^vr) is an affine fiber bundle (see Carinena et al. |t9j for the case 




Diagram 1. Chain of jets 

= 1, or Saunders [36] for the general case), which is modeled on the vector bundle 

TxU{S'T*M)®7i*{V7:), 
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where S^T*M is the space of symmetric /c-tensors on M and Vvr is the vertical fiber bundle on vr. 
Thus, taking repeated jets, {J^TCk, {T!'k)i,o, J^t^) is also an affine fiber bundle. Furthermore, J^'^^tt 
can be naturally embedded into J^vr^ (see Diagram The inclusion map ii ^ : J^^^n ^ J^Tik is 
given by 

on J^TTk (with < |/| < k), then ii,j!c( J'^^^tt) is given 



If we consider fiber coordinates 
by the equations 



a 



a 



for < |/| < A; 
when III = IJI : 



1; and 

k and / + Ij 



J + 1 



J- 



jk+i 




, '(7rfc)i,o 




M 



DIAGRAM 2. Iterated jet 



As we have said, TTk,k-i '■ J^t^ - 
dual affine bundle i^'lk-i '■ '^^''^^ ~ 
extended dual bundle {J^tx^ tt], ^ 



— y ^vr is an affine bundle, so we may consider its extended 
J^'^V and its dual affine bundle vr^ : J'^vr* — )■ J^'^ti. The 
J^~^tt) is a fiber bundle whose fibers consist of affine maps of 



the corresponding fibers of the affine bundle (J'^tt, Hk^k-i, In its turn, {J'^n*, 7^1 k-i, J^~^7i) 

is a fiber bundle whose fibers consist of classes of affine maps of the corresponding fibers of the 
affine bundle ( J^tt, VTfc.fc-i, J'^'^tt), which differ by a constant. It can be shown that there exist 
canonical isomorphisms such that J'^vr'^ ^ A™(J^^V) and J'^n* ^ K^{J^~^ii) / AJ^^J^'^-k), where 
A™(j'^~^7r) is the bundle of those m-forms over J'^-V that are annihilated when r of their 
arguments are vertical with respect to T^k-i '■ J^^^tx — > M. Locally, the elements of A^( J^~^7r) 
are of the form 

where < |/| < — 1. Thus, adapted coordinates on J'^'^tt induce coordinates on A^(J^~^7r) 
and A™(j'=-V)/A™(J^-i7r) of the form 

{x\u'^,p,pi;') and {x\u'},pi;'), 

respectively. 

While J^TT^ is naturally paired with J'^tt, A^{J^~^7c) has a canonical multisymplectic form (see 
El |9]). The pairing between the elements of J'^tt and A™(j'^~^7r) is given by 

(2.1) ^(j^(f),ujjk~i^) = a{x), such that a(a:;)?7(x) 

which is written in adapted coordinates 

(2.2) 
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where |/| = 0, . . . , k — 1 and \K\ = k. The canonical multisymplectic (m + l)-form on A™( J'^^^tt), 
which will be denoted Q, is written in coordinates 

(2.3) Q = -dpA (Tx - dpi;' A du"} A rf^-^Xi. 



3. The Skinner-Rusk formalism 

The generalization of the Skinner-Rusk formalism to higher order classical field theories will 
take place in the fibered product 

(3.1) Wo = J'7,Xj.-.^A^{j'-'7r). 

The first order case is covered in PT| [T3] : see also [3^, 130] for the original treatment by Skinner 
and Rusk. The projection on the z-th factor will be denoted pr^ (with i = 1,2) and the projection 
as fiber bundle over J^^^ti will be t^Wo,j^-"'tt = '^k,k~i ° P^i (see Diagram |3]). On Wq, adapted 
coordinate systems are of the form (x*, m", M^,p^'*,p), where |/| = 0, . . . , A; — 1 and \K\ = k. 




M 



Diagram 3. The Skinner- Rusk framework 

Assume that L : J'^tt — > M is a Lagrangian function. Together with the pairing $ (equations 
(12. ip and (12. 2p ). we use this Lagrangian L to define a dynamical function Hq (corresponding to 
the Hamiltonian) on Wq: 

(3.2) Hq = ^~Lo pr^. 

Consider the canonical multisymplectic (m + l)-form Vl on A^(j'^^^7r) (equation fl2.3p ). whose 
puUback to Wq shall be denoted also by Vt. We define on Wq the (m + l)-form 

(3.3) Qho = ^ + dHo A r]. 
In adapted coordinates 

(3.4) Ho = p'^'u^^^+p-Lix^u^^u'i,) 

(3.5) fi^o = -dp'^' A rfw? A d'^-^Xi + {p'^'du^^^ + u%^4p'^' - ^du'}^ A d'^x, 
where III = k — 1 and IJI = k. 
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3.1. The dynamical equation. We search for an Ehresmann connection F in the fiber bundle 
7TWo,M '■ Wo — > M whose horizontal projector be a solution of the dynamical equation (see 
Appendix [XI): 

(3.6) ih^Ho = {m- l)nHo- 

We will show that such a solution does not exist on the whole Wq. Thus, we need to restrict to 
the space on where such a solution exists, that is on 

Wi = {w eWo / 3h^ : TujWo — > T^Wq linear such that = h^„, 
kerh^ = {V7rwo,M)w, ih^^Hoiw) = (m - l)il/^„(w)}. 

Remark 3.1. Equation (13. 6 P is a generalization of equations that usually appear in first order 
field theories. In this particular case, from a given Lagrangian function L : J^tt — t- M we may 
construct a unique {n + l)-form (the Poincare-Cartan (n+l)-form). Hence, we have a geomet- 
rical characterization of the Euler-Lagrange equations for L as follows. Let F be an Ehresmann 
connection in tti o : J^vr — > M, with horizontal projector h. Consider the equation 

(3.8) li^nL = {n - i)nL. 

If h has locally the from 

d \ d d d 

+ Af—- + A". 



then a direct computation shows that equation (13. 8p holds if and only if 
(3-9) -^\=0 




(310) ^-^-^-A^^-^-A^^-^ + (A^-u^)-^ = 

J ^ 3 

(see [29j). If the lagrangian L is regular, then Eq. (13.91) implies that Af = uf and therefore Eq. 
(I3.10p becomes 

(3 11) ^_^_#^J^_#^!^ = 

Now, if cr(x^) = (j'^(x), af{x)) is an integral section of h we would have 

uf = — — and A - 



dx^ ^-^ dx^ ' 

which proves that Eq. (I3.1ip is nothing but the Euler-Lagrange equations for L. 

We may think Equation (13. 6p as a generalization of equation 13.81 giving the Euler-Lagrange 
equations for higher-order field theories in a univocal way, as we will see. 

In a local chart (x*, Mj,p^'*,p) of Wq, a horizontal projector h must have the expression: 

/ (9 (9 f) (9 \ 

(3.12) h = — + A" ^ + Bi)^ + a— ® dx\ 

^ ' \dx3 -^^du} ^dp/ dp J 
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where 1/1 = k — 1 and IJI = k. We then obtain that 



dL 



fc-i 



|/'|=i 



dL 



V 
fc-2 



|/|=0 



|7|=fc-l 



i+u 



dL 

\K\=k ^ 
fe-1 

|7|=0 

Equating this to zero and using Lemma IB. 31 from Appendix [Bl we have that 



Ad'^x. 



(3.13) Al 
(3.14) 



u 



0, 



Jc — 1, i = 1, . . . , m; 



5 ^. 



(3.15) 
(3.16) 



I,i 



dL 

m + 1 f dL 
/f^) + ] 



A: -2, i = l, 



dL 



1 \du 



+ Q 



li 

a I 1 



k — 1, i 



L,...,m; 



where the Q's are arbitrary functions such that 
(3.17) 



k. 



E iFTr^^^ = 0' -th|j| = i, 

I+li=J ' ' 

Remark 3.2. The ambiguity in the definition of the Legendre transform, and therefore of the 
Cartan form, becomes apparent in the equations (I3.15P and f l3.16p . as noted by Crampin and 
Saunders (see [37]). There are too many momentum variables to be related univocally with the 
velocity counterpart. To fix this, a choice of arbitrary functions Q satisfying fl3.17p must be done. 
The choice may be encoded as an additional geometric structure, like a connection. 

Applying fl3.17p to f l3.15p and f l3.16p . and using the identity ( 1B.2p . we finally obtain the equa- 



tions 




(3.18) 




(3.19) 




(3.20) 


E Pi" 




I+U=J 


(3.21) 


E Pi' 



with |/| = 0, 

dL 



dW^' 
dL 



I+li=K 



du'} 

dL 
dv^' 



B^, with \J\ 



k — 1, i 



l,...,k-l; 



with li^l = k. 



Notice that equation (13.2ip is the constraint that defines the space Wi, and that (13.18p . (13.19P 
an d fl3.20p are conditions on coefficients of the horizontal projectors h. Note also that, for the 
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time being, the A's with greatest order index and the C"s remain undetermined, as well as the 
most part of the i?'s. From the definition of Wi, we know that for each point w G Wi there 
exists a horizontal projector : T^Wq — > T^Wq satisfying equation (13. 6p . However, we cannot 
ensure that such h^, for each w G Wi, will take values in T^Wi. Therefore, we impose the natural 
regularizing condition h^{T^WQ) C T^Wi, Vw G Wi. This latter condition is equivalent to having 



PL 



0, 



which in turn is equivalent (using (I3.12p and (13.181) ) to 



with \K\ = k. Thus, if the matrix of second order partial derivatives of L with respect to the 
"velocities" of highest order 



(3.23) 



dujdu'j^ 



is non-degenerate, then the highest order A's are completely determined in terms of the highest 
order S's. In the sequel, we will say that the Lagrangian L : J^tt — > M is regular if, for any 
system of adapted coordinates the matrix, (I3.23P is non-degenerate. 

Up to now, no meaning has been assigned to the coordinate p. Consider the submanifold W2 of 
Wi defined by the restriction Hq = 0. In other words, W2 is locally characterized by the equation 

As before, we cannot ensure that a solution h of the dynamical equation ( 13. 6 p takes values in TW2- 
We thus impose to h the regularizing condition h^(T^PFo) C TwW2, Vw G W2, or equivalently 
h.{d/dx^){Ho) = 0. Therefore, the coefficients of the linear mapping h are governed by the 
equations IKWf . ( IXTOj) . (K22\f and in addition 

(3-24) ^^ = 0^ + ^%d^^ - ^^^^^^ - ^-^^^V 

Note that, thanks to the Lemma [B.3I and equation (I3.2ip . the terms with A's with multi-index 
of length k cancel out, and the A's with lower multi-index are already determined. So, in some 
sense, the C"s depend only on the 5's. 

3.2. Description of the solutions. The relations (l3.2Up (with | J| = k — V) and (I3.22p can be 

seen as a system of linear equations with respect to the i?'s. When /c = 1, equation (I3.19P should 
be considered instead of equation (I3.20p . In the following, we are going to suppose that n = 1, 
since the dimension of the fibres is irrelevant for our purposes and we may ignore it. The number 
of 5's with order /c — 1 (with multi-index length k — V) is given by 



m — 1 + k — 1 
m — 1 
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and the number of equations with such 5's is 

m — 1 + k\ /m — l + Zc — 

■ m + 

m — I J \ m — I 

An easy computation shows that the system is overdetermined if and only if /c = 1 or m = 1 
(examples 14.11 and I4.2p . and completely determined when k = m = 2. In all other cases the 
system is underdetermined, but it still has maximal rank. 

Proposition 3.3. Suppose that k > 2 and m > 2. Then, the system of linear equations with 
respect to the B 's 

(3.25) Ei'? - Z/'N 

I+li=K " |/|=0 \J\=k 

where | J| = A; — 1, j = 1, . . . , m and \K\ = k, has maximal rank. 

Proof. In a first step, we are going to describe how to write the matrix of coefficients. Then, we 
will select the proper columns of this matrix to obtain a new square matrix of maximal size. We 
finally shall prove that this matrix has maximal rank. 

The matrix of coefficients will be a rectangular matrix formed by I's and O's. The columns will 
be indexed by the indexes of the -B's, and the rows by the indexes of the first partial derivatives 
that appear in the equations fl3.25p and f l3.26p . As B^^- has three indexes, the columns of the 
matrix of coefficients will organized in a superior level by the index i, in a middle level by the 
index j and in an inferior level by the multi-index I. The rows will be organized at the top by 
the index J for the first equation, fl3.25p . and at the bottom by the index j and then by the 
multi-index K for the second equation, f l3.26p . 

As the matrix of coefficients has more columns than rows, we shall build a second matrix that 
has as many columns and rows as the matrix of coefficients has rows. To do that, we select a 
column of the matrix of coefficients for each row index using the following algorithm (for the sake 
of simplicity): 

01 ForEach (j,K) 

02 Define G={(I,i) :I+l_i=K} 

03 If Cardinal (G)=l 

04 Select the column (i,j,l) 

05 Elself 

06 Select a column (i,j,I) such that (I,i) is in G and i\neq j 

07 Endlf 

08 EndFor 

09 ForEach J 

10 If J(l)=k-1 

11 Select the column (m,m,J) 

12 Elself 

13 Select the column (1,1, J) 

14 Endlf 
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15 EndFor 

Now, this matrix being defined and since it is full of O's and has only few I's, we are going to 
develop its determinant by rows and columns. Notice that the columns selected at line 6 have 
only one 1 each, thus we can cross out the rows an columns related to these I's. Now the rows at 
the bottom part of the remaining matrix (related to the second equations) have only one 1 each, 
thus we can also cross out the rows an columns related to these I's. Now, the remaining matrix 
has the property of having only one 1 per column and row (there must be at least one 1 per row 
and column, and no two I's may be at the same row or column), thus its determinant is not zero 
and the matrix of coefficients has maximal rank. □ 

Another way to interpret the tangency condition (13.221) is the following one: Let us suppose 
we are dealing with a first order Lagrangian (example 14. ![ equation f l4.8p ). One could apply the 
theory of connections to the Lagrangian setting and the Hamiltonian one as separate frameworks. 
We know that they must be related by means of the Legendre transform and so are the horizontal 
projectors induced by these connections. Thus, equation (14. 8p is nothing else than the relation 
between the coefficients of these horizontal projectors. 

3.3. The reduced mixed space W2. In section §3.21 we reduced the space Wi to W2 by consid- 
ering the constraint Hq = 0, which is a way of interpreting the coordinate p as the Hamiltonian 
function. But W2 is not a mere instrument to get rid off the coordinate p or the coefficients 
Cj. As the premultisymplectic form ^Ihq, it encodes the dynamics of the system and, when L is 
regular, it is a multisymplectic space. Indeed, when k = 1, W2 is diffeomorphic to J^tt (c/. de 
Leon et al. [31j), which is not true for higher order cases. 

Proposition 3.4. Let W2 = {w G Wi : Hq{w) = 0} and define the (m + l)-form, Q2 o-s the 
pullback ofQfjo ^2 by the natural inclusion i : W2 Wq, that is Q2 = '^*(^//o)- Suppose that 
dimM > 1, then, the (m + l)-form Q2 is multisymplectic if and only if L is regular. 

Proof. First of all, let us make some considerations. By definition, ^2 is multisymplectic whenever 
has trivial kernel, that is. 




This is equivalent to say that 




Let V G TWo be a tangent vector whose coefficients in an adapted basis are given by 



_|_ /1Q _|_ _|_ 

dx' ^du^j ° dpi: dp 



Using the expression (13. 5p . we may compute the contraction of by v, 

i^riHo = -Bl'duJ A d'^-^Xi + AJdp^^ A d'^-^Xi - X^dpi' A duj A d'^-^x, 
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On the other hand, if we now suppose that v is tangent to W2 in Wo, that is f G i^{TW2), we 
then have that 



(3.28) 



d 



( E 

\I+U=K 



dL 



du 



{v) = and dHo{v) = 0, 



K 



which leads us to the following relations for the coefficients of v, 



(3.29) 
(3.30) 



I+li=K 



dx^du 



K 



^ dujdu%- 



and 



Aa li I Uli^.a 



\i. dL /[a dL 



0. 



It is important to note that thanks to Lemma IB. II and the equation fl3.2ip which defines Wi (and 
hence W2)-, the terms in (13.271) and (I3.30p involving A's with multi-index of length k cancel each 
other out. 

These considerations being made, suppose that ^2 is multisymplectic and, by reductio ad 
absurdum, suppose in addition that L is not regular, which means that the matrix 



has non-trivial kernel. Let v G TWq be a tangent vector such that all its coefficients are null 
except the A's of highest order which are in such a way they are mapped to zero by the "hessian" 
of L. Such a vector v fulfills the restrictions fl3.29p and fl3.30p . thus it must be tangent to W2 in 
Wq, V G i^{TW2)- But, as iyflno has no A's of highest order, it must be zero, i^Qno = 0; which 
is a contradiction. 

Conversely, let us suppose that L is regular, then equation (13.21 p defines implicitly the coordi- 
nates as functions of the other coordinates. That is, locally there exist functions /^(a;*, m", p^'*) 
such that = on i{W2). Furthermore, given a system of adapted coordinates (a;*, m", m^, p^'*, p) 
on Wq, (a;*, Mj,p^'*) defines a coordinate system on Wq and the inclusion is given by: 

k-2 



1/1=0 



|/|=fc-i 



From equation (13. 5p . we can compute an explicit expression of the (m -|- l)-form in this 
coordinate system. 



fc-i 



I\=0 
k-2 



+ 



+ 



k-l 



dL 



|/|=0 1/1=0 



Ad"'x 



E {p'»'dfr^u + f?+iM') 
\i\=k-i 



\K\=kI+li=K 



Ad'^X, 
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where we have used equation (I3.2ip in the last term. Note that, by Lemma [B.lt the first and last 
terms of the last bracket cancel each other out. Now, 

ig/Q^^n^ = dpi;' A du? A d"'-'^Xij -[...] A d^'-^xj 



--1 J 



IqIq^iM^ = du" A d'^ 'xi + 

where < |/| < A; — 1. We deduce from here that the kernel of Q2 is trivial, kerf22 = {0}, and 
is multisymplectic. □ 

Note 3.5. In the particular case when dimM = 1, the Lagrangian function L : J'^vr — > M is 
regular if and only if the pair {fl2,Tw2 udt) is a cosymplectic structure on W2- We recall that a 
cosymplectic structure on a manifold of odd dimension 2n + 1 is a pair which consists of a 
closed 2-form Q and a closed 1-form 77 such that 77 A fi" is a volume form. 

We remark that, if the Lagrangian L is regular or (from Proposition 13. 3p if A;, m > 1, then there 
locally exist solutions h of the dynamical equations f l3.6p on W2 that give rise to connections F 
in the fibration ttw^m '■ — > M along the submanifold W2 (see Appendix [A]). In such 
global solution is obtained using partitions of the unity, and we obtain by restriction a connection 
r, with horizontal projector h, in the fibre bundle vtvi/jM '■ W2 — M, which is a solution of 
equation (13. 6p when it is restricted to W2 (in fact, we have a family of such solutions). 

In some cases, but only when dimM = 1 or A; = 1, it would be necessary to consider a subset 
defined in order to satisfy the tangency conditions fl3.22p and fl3.24p : 



W 1 



W3 = {w eW2 / 3h^ -.T^Wo — > T^W2 linear such that = h 
kerh^ = {Vttwom)'^^ ih^^Hoii^) = {m - l)fi^o(u;)}. 

We will assume that is a submanifold of 1^2- If h^(T^Wo) is not contained in T^Ws, we go 
to the third step, and so on. At the end, and if the system has solutions, we will find a final 
constraint submanifold Wf, fibered over M (or over some open subset of M) and a connection 
F/ in this fibration such that Fj is a solution of equation (13. 6 P restricted to Wf. 

In any case, one obtains the Euler-Lagrange equations. In the following result, Wf denotes 
the final constraint manifold, which is W2 when k,m > 1, and h the horizontal projector of a 
connection in 7rw2,M '■ Wf — > M along Wf, which is solution of the dynamical equation. 

Proposition 3.6. Let a be a section of TTWf^M '■ Wf — > M and denote a = i o a, where 
i : Wf Wq is the canonical inclusion. If a is an integral section of h, then a is holonomic, in 
the sense that 

(3.31) prioa = f{TTwj,E'=>^), 
and satisfies the higher-order Euler-Lagrange equations: 

(3.32) ,-K,^,„.).(t,_l)U,^|£).0. 



|J|=0 
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Proof. U a = (x*, cTj, o"^'*, a) is an integral section of h, then 



dx^ dx^ ™ dx^ 



where the A's, S's and C"s are the coefficients given in (I3.12p . From equation (I3.18p . we have 
that cr is holonomic, in the sense that crf^^ = daj jdx^ . On the other hand, using the equations 
( I3.19p . (I3.20p and f l3.2ip . we obtain the relations (where = pri o a): 

(3.34) ^^'^ = 9^°^^- S^T' -i^h|^l = l'---'^-l' 
/+ii=J 

dL 

(3.35) 5^ a^'^ = J— o (j), with \K\ = k. 



I+U=K 



K 



From the equations (I3.33P and f l3.34p for \J\ = 1 we get 

aw" ox^ 



\I\=1 ^ ^ |/|=1 I 



Applying now Lemma [B. II on the last term and repeating this process until |/| = — 1 we reach 



- ovr|^-Eo>r(|;||)^E|^E-? 

I ri_i \ J ^ I Ti—o r_i_i .— 7 



|/|=1 ^ ^ |J|=2 I+li=J 



Qr^a ^ (dx^du'f) dx^ dx' 

\I\=1 ^ ^ |/|=2 ^ ^ |/|=2 i 



|/|=0 ^ ^ ^ |/|=fe-l i 

EM)"'o"'^r(^l9-(-i)'-'E|:;^ E 

1/1=0 ^ ^ |i^|=fc I+U=K 



where by abuse of notation = j^^\'KWf,E ° cr)- Finally, it only rest to use equation fl3.35p to 
prove the desired result. □ 



4. Examples 

First, we are going to study the particular cases when k = 1 and m = 1, which correspond to 
the First Order Classical Field Theory and to the Higher Order Mechanical Systems, respectively. 
Theoretic results for these cases are very well known [H [3ll [131 [23] ^^'^ we are only going to recover 
these results from our general setting. 
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Example 4.1 (First order Lagrangians {k = 1)). Let us suppose that k = 1, which corresponds 
to the case of first order Lagrangians. In that case the velocity-momentum space is Wq = 
J^TT with adapted coordinates {x^,u'^,u1,p,pl^). The premultisymplectic (m + l)-form 

would be 

(dL 8T \ 

pidut + ufdpi - ^du'' - ^dutj A (Tx, 

and horizontal projectors on TWq would have locally the form: 

f d d d d d\ 

(4.2) h = — + + A%—- + B^— + Cj— ® dx'. 

^ ' \dx'J ^ dW^ '^duf ^dpl^ ^dpj 

Solutions of the dynamical equation would satisfy the relations 



(4.4) = for z = 1, . . . ,m; 

(4.5) At = <, for2 = l,...,m; 
from which we deduce the Euler-Lagrange equations 

(4.6) f{Trw2,M°(^) 
where W2 is defined by 




(4.7) W2 = \ {x\u-,ut,p,pl) eW, : pl = -^, p = L- Y,P'u, 



dut' 

' 1=1 



We then obtain the tangency conditions: 



(48) B^ - -^ + u^^^ + TA^ ^'"^ 



I — 1 t ^ 



(4-9) = ^ + 



Note that (14. 8 P is the relation that would appear between the coefficients of a Lagrangian and 
a Hamiltonian setting through the Legendre transform. For simplicity, suppose that n = 1 and 
ignore the a's and /3's that appear above. Consider the linear system of equations with respect 
to the B's formed by equations (14. 3 p and (14. 8p . This system is overdetermined since it has -|- 1 
equations and only variables [Bj). 

Example 4.2 (Higher order mechanical systems (m = 1)). Let us suppose that m = 1, which 
corresponds to the case of mechanical systems. In that case the velocity-momentum space is 
Wo = J^n X jfe-ijr A™(j'^~^7r). Since here a multi-index J is of the form (/) with 1 < / < /c, we 
change the usual notation for coordinates to 

u^j^ufj^ and p^'i— 
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and we adapt the remaining objects to this notation. So adapted coordinates on Wq are of the 
form {x, u'^,p,Pa), where I = 1, . . . , k. The premultisymplectic (m + l)-form would be 

(4.10) Q^^ = -Y^ dp'+' Adu^ + J2 {Pad< + Adx-Y, g^du^ A dx, 

1=0 1=1 1=0 ' 



and horizontal projectors on TWq would have locally the form: 

Solutions of the dynamical equation would satisfy the relations 
(4.12) Bi 



d_ 

dp 



dx. 



(4.13) = in^-B'^\ for/ = l,...,A;-l; 



(4.14) p 



k 



duf 
dL 



(4.15) = for / = 0,...,A;- 1. 

which we use to get the Euler-Lagrange equations 

(4.16) 

\i=o ' 



where W2 is defined by 



dL 



(4.17) W2={{x\u",ur,p,pl)eW, : = p = L -J^Pa^^T 



, p = 1^- 2_^ f^ 

"'^ 1=1 

We then obtain the tangency conditions: 



Q2r ^=-1 par Q2r 

(4.18) Bl = + y ul^^U^ + Ai^^ = 

^^^^^ ^ du^.dul ^ duidul 

1=0 ' ^ 1=1 

Note that, thanks to equation (I4.14p . the terms in (I4.19P with coefficient A^ cancel out. Now, 
for simplicity, suppose that n = 1 and ignore the a's and /3's that appear above. Consider the 
linear system of equations with respect to the i?'s formed by equations fl4.13p (with I = k — 1) 
and f l4.18p . This system is overdetermined since it has 2 equations and only one variable (B^). 

Example 4.3 (The loaded and clamped plate). Let us set M = and E = x M = M^, and 
consider the Lagrangian 

L{X, y, U, Ucc, Uy, Ur,^, U^y, Uyy) = ^{ul^ + 1u\y + " 2^^), 
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where q = q{x, y) is the normal load on the plate. Given a regular region R of the plane, we look 
for the extremizers of the functional I{u) = J^L such that u = du/dn = on the border dR, 
where n is the normal exterior vector. The Euler-Lagange equation associated to the problem is 

(4.20) 

Written in the multi-index notation, the Lagrangian has the form 

L{f4>) = ^(^(2,0) + 2m?i,i) + M(o,2) - 2gu) 
and the Euler-Lagrange equation reads 

M(4,0) + 2M(2,2) + M(0,4) = Q- 

The velocity-momentum space is Wq = J^tt Xji^ A2(J^7r), with adapted coordinates {x,y,Ux, 
Uy,Uxx,Uxy,Uyy,p,p^,p'^,p^y,p^y,p'^^,py^). It Is st r alght f or War d to write down the premultisym- 
plectic 3-form and a general horizontal projector on TWq, so we are not going to do it here. Even 
so, the coefficients of solutions of the dynamical equation would satisfy the relations 

/y-.XX fi , 

(4.21) Bf + B/ = -2q J Z j{y,xl Jy,y p-v + pV- = 2u^y 

P JD^ JDy yy _ 

P Uyy 

where the latter ones are the equations that define Wi. The tangency condition on Wi gives us 
the relations 

T^XjX A A 

J->x — ■'^xx,x -'-'y ^xx,y 

(4.22) Bl^y + Bl'^ = 2Axy^x B^^v + Bf^ = 2A.,,, 

Dy,y ^ A r>y,y ^ A 

^x ^yy,x ^y ^yy,y 

from where we can see that the Lagrangian is "regular", since 



(4.23) 



\K\ = \K'\=2 




Finally, we remark that the middle equations of (14.211) and (I4.22p form a 8 x 8 linear system of 
equations on the i?'s, which is completely determined. 

Example 4.4 (The Camassa-Holm equation). In 1993, Camassa and Holm introduced the following 
completely integrable bi-Hamiltonian equation (see [5j): 

(4.24) Vt - Vyyt = -3VVy + 2VyVyy + VVyyy , 

which is used to model the breaking waves in shallow waters as the Korteweg-de Vries equation. 
But, as the former is of higher order, we are going to use it as example. 

The CH equation fl4.24p is expressed in terms of the Eulerian or spatial velocity field u{y,t), 
and it is the Euler-Poincare equation of the reduced Lagrangian 

(4.25) /(t;) = ^ j {v' + vl)dy. 

To give a multisymplectic approach to the problem, as Kouranbaeva and ShkoUer did (see |26]). 
we must express the CH equation (14.241) in Lagrangian terms. Thus, we shall use the Lagrangian 



18 



C. M. CAMPOS, M. DE LEON, D. MARTIN DE DIEGO, AND J. VANKERSCHAVER 



variable u{x,t) that arises as the solution of 

(4.26) ^^ = v{u{x,t),t). 

The independent variables [x, t) are coordinates for the base space M = S*^ xR, and the dependent 
variable u{x,t) is a fiber coordinate for the total space ii^ = S*^ x M x M = S*^ x M^. The Lagrangian 
action is now written as 

(4.27) L(x, t, u, Ux,ut, u^x, Uxt, uu) = ^ {u^uf + u~'^ult) 

The coefficients of a horizontal projector which is solution of the dynamical equation must satisfy 

B^'^^ + B/ = 

= l/2{ul-{uxt/uxf)-{B:^^ + Bf) 

pxx ^ g 
ptt ^ Q 

where the last three are the equations that define Wi. The tangency condition on Wi gives us 
the relations 

^X,X ^ Q 

^x,t _^ Qt,x ^ -U^^^^UxxUxt + Axt,xU~^ 
J^t,t g 

(4.29) j^f,x I Q 

^x,t _^ ^t,x ^ _(y^^/^^)2 _^ Axt,tU~'^ 

bY = 

from where we can see that the Lagrangian is clearly "singular", since 

/ <92/ \ /O 

(4.30) ={ 

\OUkOUk' J \K\ = \K'\=2 \ 

Again, we may form a completely determined system of linear equations on the i?'s with the 
corresponding relations of fl4.2ip and the equations (I4.29p . 

Example 4.5 (First order Lagrangian as second order). For the sake of simplicity, let suppose that 
n = 1. Given a first order Lagrangian L : J^tt — > M, extend it to a second order Lagrangian, L = 
L o 7r2,i. Consider the first and second order velocity- momenta mixed spaces Wq = J^vr x e A'^E 
and Wq = J^vr x A^( J^vr), with adapted coordinates {x\u,Ui,p,p^) and {x\u,Ui,UK,P,p\p^-') 
(with \K\ = 2), respectively. Let tTq'^ : — > Wq be the natural projection (DiagramHj). 

We are going to apply the theory we have developed here to the systems given by each La- 
grangian. Consider the premultisymplectic forms Qhq ^i-nd ^h^, where Hq and Hq are the corre- 
sponding dynamical functions (equations fl3.2p and (13. 3p ). Let h and h denote solutions of the 
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Diagram 4. The 1st and 2nd order Lagrangian settings 

respective dynamical equations on {Wq,Qho) {WQ,Qffg). They would locally have the form 

d d d ■ d d \ 

+ Ai— + An— + Bl— + d— O dx\ 



dx^ ^ du dui ^ dp^ ^ dp 



(4.31) B] 

(4.32) p' 



^ dx^ du dui duK dp^ dp^"^ dp 

where \K\ =2. We then obtain the relations 

du 
dL 

dui 

(4.33) A = Ui, 
for {W^,nHo,'^); and 

dr 

(4.34) - ^, 

dr 

dr 

(4.36) p'^+p^' = (l. + l,)!-T^ = 0, 



+1, 



(4.37) Ai = Ui, 

(4.38) Aij = ui^+i., 

for {W^,nH^,h). Equations and f lOB]) . together with i^o = and i?o = 0, define the 

corresponding submanifolds W2 and W2 of W^q^ and Wq. 

We notice that, even though L is in some sense the same Lagrangian than L, a solution of 
the dynamical equation on Wq may be easily determined, while in Wq the space of solutions has 
grown (there are more coefficients to be determined). We thus infer from here, that a solution 
h of the dynamical equation in Wq must satisfy an extra condition. Since p = L — p^Ui + in 
W2, the projection tTq'^ maps W2 to Wj^. We therefore impose to a solution h of the dynamical 
equation along W2 to be in addition projectable to a solution h of the dynamical equation along 
W}. In such a case, we would have that 



(4.39) B 
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which imphes that the following equation 

rlT 

(4.40) ,f = ^ 

is now a restriction in So, by tangency condition, we get 

^ dx^dui dudui ^^''^^^ dukdui dx^ dut 

Combining this with equation (I4.3ip . we finally obtain 

f4 42) d dL ^ 

du dx^ duj 

which is the Euler-Lagrange equation. 

It is worth to remark here that, at this time, the Euler-Lagrange equation has not been deduced 
by the process shown in the proof of Proposition [XHl but directly from the projectability condition, 
although the previous Euler-Lagrange equation may be recovered from any of the two settings. 

5. Conclusion 

We have developed an intrinsic and global expression for the Euler-Lagrange equations for 
higher-order field theories. The main ingredients of this setting are the mixed space of velocities 
and momenta Wo = J^tt Xjk-i^ A™(j'^^^7r) and the premultisymplectic form 

^Iho = ^ + dHo A T] 

defined on it, which encodes the dynamics of the system through the dynamical equation 

We have analyzed in detail the existence of solution of this equation. Our approach gives rise to 
an unambiguous formulation of Lagrangian field theories of higher order. 

In a future paper we will explore the extension of our techniques to the case of higher-order 
field theories with constraints, optimal control problems for partial differential equations and the 
implementation of numerical methods obtained directly from our approximation. 

Appendix A. Connections 

A connection F in a fibration t^p^m '■ P — > M is given by a horizontal distribution H which is 
complementary to the vertical one Vttp^m, that is 

TP = H © Vttp^m- 

Associated to the connection there exists a horizontal projector h : TP — y H defined in the 
obvious manner. If (x',?/") are fibered coordinates on P, then H is locally spanned by the local 
vector fields 

d V d .a 



^dx'- J * ' dy"- ' 

{d/dx^)^ is called the horizontal lift of d/dx^, and are the Christoffel components of the 
connection. 

Assume that np'^M '■ P' — ^ ^ and t^p^m '■ P — ^ ^ are two fibrations with the same base 
manifold M, and that T : P' — t- P is a surjective submersion (in other words, a fibration as 
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well) preserving the fibrations, say, t^p^m o T = iTp/^M (Diagram |5]). Let T' be a connection in 
TTpi^M '■ P' — > M with horizontal projector h. 

T 

P' 




M 

Diagram 5. Preserved fibration 

Definition A.l. V is said to be projectable if the subspaces TT{z'){H.z') do not depend on 
z' e T-\T{z')). 

If r' is projectable, then we define a connection F in the fibration ttp^m '■ P — ^ ^ as follows: 
The horizontal subspace at z € P is given by 

=TT(^')(H.O, 

for an arbitrary z' in the fibre of T over z. It is routine to prove that H defines a horizontal 
distribution in the fibration ttp^m '■ P — > M. 

We can choose fibered coordinates {x'^,y"',z") on P' such that {x\y"') are fibered coordinates 
on P. The Christoffel components of F' are obtained by computing the horizontal lift 

f d \^ d d d 

(a?j =S;r + r?(..y..)^ + rr(..y..)— . 

A simple computation shows that F' is projectable if and only if the Christoffel components F" 
are constant along the fibres of T, say F" = F?(a;, y). In this case, the horizontal lift of djdx^ 
with respect to F is just 

As an exercise, the reader can easily check that, conversely, given a connection F in the fibration 
T^PM '■ P — ^ M and a surjective submersion T : P' — > P preserving the fibrations, one can 
construct a connection F' in the fibration ttp'^m '■ P' — > M which projects onto F. 

The notion of connection in a fibration admits a useful generalization to submanifolds of the 
total space. Let ttpm : P — > M be a fibration and N a submanifold of P. 

Definition A. 2. A connection in ttp^m '■ P — > M along the submanifold N consists of a family 
of linear mappings 

K : T,P T,N 
for all z & N, satisfying the following properties 

hl = hz, kerhz = V71p^m)z, 

for all z & N. The connection is said to be differentiable (flat) if the distribution imh C TN is 
smooth (integrable). 

We have the following. 

Proposition A. 3. Let h be a connection in ttp^m '■ P — > M along a submanifold N of P. Then: 
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(1) 7Cp^m{N) is an open subset of M. 

(2) (vrp,M)|Ar : N — > np^M{N) is a fibration. 

(3) The 1 -jet prolongation J^(7rp^j\/)|^ is a submanifold of J^tip^m- 

(4) There exists an induced true connection Fat in the fibration {'n'pM)\M '■ N — > ttpm{N) 
with the same horizontal subspaces. 

(5) is fiat if and only ifh is fiat. 



Proof See [301131]. □ 

Appendix B. Multi-index properties 

This section is devoted to some simple, but useful, properties of multi-indexes. 

Lemma B.l. Let {aj^i}j ^ be a family of real numbers indexed by a multi-index I G N"* and by 
an integer i such that 1 < i < m. Given an integer I > 1, we have that 

m 

(B.l) E E«^.^ = E E 

\i\=i-ii=i \j\=ii+ii=j 

Proof. The proof is trivial when we realize that the sets {(/,«) : |/| = / — 1, 1 < i < to] and 
{(J,2) : / + 1, = J, I J| = /} are in bijective correspondence. □ 

Lemma B.2. Let J G N™ be a fixed multi-index. We have that 

i+ii=j I I ^ 

Proof. 

^ \J\ ^ \A ^ 1^1 + 1 

□ 

Lemma B.3. Let |aj, &"^} j be a family of real numbers indexed by a multi-index J G N™. Given 
an integer / > 1, we have that 

(B.3) E^'«^= E E^^(^'^''+^''>^+i^' 

\j\=i i/i=/-i «=i 

where {Q'^'*}^ ^ family of real numbers such that for any multi-index J G |J| > 

we have that 

I+li=J I I 

Proo/. 



\j\=ii+ii=j 
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